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Abstract 

'&) , 

We consider three methods by which one can generate new cosmological models. Two of these are 

o' 

£NJ ' based on the Lorentzian structure of spacetime. In a Lorentzian manifold there can exist horizons 

that separate regions of spacetime that can be interpreted as cosmological models from others that 
have the character of "black holes." A number of well known solutions of this type can be used to 
generate both known cosmological models and others that do not seem to have been recognized. 
Another method based on the Lorentzian character of spacetime is to simply interchange some 
space variable with time and try to restructure the metric to make a viable cosmology. 

o 

A more broad-ranging method is the use of modern solution-generating techniques to construct 

new models. This method has been widely used to generate black hole solutions, but seems not to 

O ■ have been so widely used in cosmology. We will discuss examples of all three methods. 
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I. INTRODUCTION 



The standard model of cosmology assumes that the universe today is homogeneous and 
isotropic, which means that there are six Killing vectors of the manifold that give the 
isometries that realize these symmetries. However, in earlier stages of the universe it is 
assumed that there might have been large amounts of anisotropy and inhomogeneity that, 
by some physical process, could have been reduced to the point where today we observe 
the standard model (for a compendium of reasons for considering inhomogeneous models, 
see Ref. lj). Anisotropic and inhomogeneous metrics have fewer and fewer Killing vectors, 
and we could, in principle, arrive at a completely general solution of the Einstein equations 
as a cosmological model. At some point it will become difficult to distinguish between a 
cosmological model and any other general solution of the Einstein equations, except as a 
matter of interpretation. In fact, even an inhomogeneous cosmology with two spacelike 
Killing vectors is difficult to distinguish (locally) from gravitational waves propagating in 
one space direction. One method that has been used to separate cosmological models from 
pure wave solutions is the prescription of their global topology. If we insist that cosmologies 
have compact t = constant surfaces, then a number of gravitational wave solutions can be 
made into cosmologies by compactifying in certain directions of the system. We will give 
examples of this procedure below. If we allow open universes, whether we call a solution a 
cosmological model or not is again a matter of interpretation. 

With these caveats we can give a number of methods by means of which one can transform 
known solutions of the Einstein equations into cosmological models. These methods can be 
broken down into three broad classes: 

1) Horizon methods 

2) Causal structure methods 

3) Mapping methods 

The first of these methods has a long history, even though it never seems to have been 
thought of as a "method." Perhaps the oldest example is the deSitter metric Q|, although it 
is a somewhat degenerate example. Originally deSitter proposed a solution to the Einstein 



equations with cosmological constant A of the form 



3, 
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ds> = - (l - ±f\ dt> + *^ + r\df + sin 2 6d^). (1) 
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This seems to be a static metric similar in form to the Schwarzschild metric, and even has 



a "singularity" (which caused much comment at the time) at r = J 3/ A. Of course, this 



singularity is just a horizon where the light cones tip over sufficiently that for r > y 3/A, 
d/dr becomes timelike and d/dt spacelike. For large r, then, we can rename the coordinates 
(r — ► t, t — *■ f), and the new metric is 

ds 2 = -^|r— j" + " X ) df2 + + sin2 ( 2 ) 

which is an obvious cosmological model. It was "recognized" (the reason for using the quota- 
tion marks will become obvious below) relatively quickly that this was indeed a cosmological 
model by making use of the coordinate transformation 



t-rW**, r = T --^-^ 1 

v sin 6 cos (p = X, r sin 6 sin = Y, rcos9 = Z, (3) 

which makes the metric 

ds 2 = _ dT 2 + e 2^T^ x 2 + rfy2 + ^ 2]) (4) 

which is a = Friedmann-Robertson- Walker metric. 

Here we should study the Penrose diagram of this metric. In Figure ^ we give a Penrose 
diagram that will cover several cases that we will discuss later. Because of this we will give 
only generic values in the figure, and each case will correspond to different values of these 



parameters. In the deSitter case, the lines aa and bb correspond to r = J 3/A, and both ±7 to 
r = 0. In region I the spacelike Killing vector d/dr generates spacelike surfaces that are the 
t = constant surfaces of (UJ. in region II the same Killing vector is now d/dt (timelike) and 
generates timelike hypersurfaces r = constant. The horizon splits the region comprised by 
both regions I and II in two. Region I might be called a "black hole" region since the metric 
is reminiscent of the Schwarzschild metric, and has the spacelike Killing vector d/dr, while 
region II is a "cosmological" region with metric (jlj). In the cosmological region there seem 
to be "singularities" at T = ±oo, but we can see that [note that T = f + \\f\ In — lj] 
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FIG. 1: A generic Penrose diagram which will be used for several metrics of the article 

the singularity at T = — oo is just the horizon at t — y3/A, so it is nothing more than a 
null surface separating two parts of the spacetime. 

The problem here is that the deSitter metric has such a high degree of symmetry (a 
space of constant curvature with 10 Killing vectors) that the distinction between regions I 
and II is artificial. In fact, the coordinate transformation (jHJ) (with t — > r and f — ► t and 
4^ — 1 — > 1- g r ' 2 ) is equally valid in region I, and, in fact, this was what was recognized by 
Lemaitre Q] and Robertson in transforming the deSitter metric into the form that was 
eventually used in the steady-state model. Actually, even Minkowski space (also a space of 
constant curvature with 10 Killing vectors) can be looked at in the same way, as we will see 
in Sec. 2. 

In Sec. 2 we will consider several metrics where horizon methods can be used to construct 
cosmological models that do not have the disadvantages of the deSitter metric, that is, there 
exists no coordinate transformation that is valid for the static region which transforms it 
into a cosmological model. The fact that "singularities" in the new cosmological models 
may just be null surfaces where one passes from one coordinate patch to another, will, of 
course, still be a feature, although true curvature singularities may also exist. One point 
that should be mentioned is that most of the models generated by this method that we will 
consider are vacuum models. There are some electrovac models that can be thought of as 
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being generated by horizon methods that we will mention, but it seems not to be known 
which matter-filled cosmological models can be generated by this method. 

The second class of methods, causal structure methods, is similar to the previous one, 
but there is no horizon to distinguish different regions of spacetime. The paradigm of this 
method is the unpolarized Gowdy T 3 vacuum spacetime p which will be discussed in more 
detail below. The idea behind this method is to begin with a known solution and simply 
rename some spacelike variable as t, and the timelike variable as a spacelike variable (we 
can call this t «-» r), and then change the causal structure of the resulting metric so that 
the a vector tangent to the "time" direction is truly timelike, and one tangent to the new 
"space" direction is spacelike. There are a number of problems with this approach. The 
first is that determining the causal structure change is not necessarily trivial. In a diagonal 
metric this difficulty can usually be handled by simply changing the signs of two of the 
metric components, but in a more general metric the procedure may be more complicated. 
The most important difficulty, however, is that the new metric is not guaranteed to be a 
solution to the Einstein equations. For matter filled models, where physical quantities can 
be used to define spacelike and timelike surfaces, the new metric will almost certainly not 
be a solution. For vacuum metrics it is difficult to tell. As with horizon methods, the known 
solutions of this type are vacuum wave solutions, and they have rather simple structures, 
and give inhomogeneous cosmological models. Here, of course, one has the problem of 
distinguishing between gravitational wave solutions and inhomogeneous models. 

A very simple example of this procedure is the plane wave |2( 

ds 2 = L 2 {u){e 2 ^dx 2 + e- 2 ^dy 2 ) + dz 2 - dt 2 , (5) 

where u = t — z. If we make the coordinate change t — > z, z — > t and change the signs of 
dz 2 and dt 2 , we find (u = t — z) 

ds 2 = L 2 (-u)[e 2 ^dx 2 + e~ 2 ^dy 2 } - dt 2 + dz 2 . (6) 

Since L and j3 are arbitrary function of their arguments, we can remove the minus signs 
in the arguments of L and (3 in (JBJ) and obtain (j3J) again. Because of the high symmetry 
of the problem, this procedure gives the same metric and Einstein equations, and there 
is no question that any solution of the Einstein equations for (j3J) gives a solution to ©. 
Here, of course, one runs into the problem of defining a cosmological model. Is either of (JHJ) 
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or (jHJ) a gravitational wave or an inhomogeneous cosmological model? Either one is both, 
depending on how we interpret them. If, however, we insist that a cosmological model have 
compact t = const, sections, then we can compactify (JHJ) in the space directions by just 
making < x, y, z < In and identifying the points at and 2ir, giving the manifold a T 3 
topology. This changes the boundary conditions on the functions L and (3 in that one must 
have L(t, 0) = L(t, 2ir) and /3(t, 0) = (3(t, 2n). Since the only Einstein equation is 

L"(u) + [P'{u) 2 ]L = 0, (7) 

' = d/du, (3 is arbitrary except for the boundary conditions, and L obeys a simple equation 
that has solutions that satisfy the boundary conditions. Here one must be careful not to 
change coordinates back to the original metric, since we will generate a spacetime with 
closed timelike lines, a problem with cosmological models generated in this manner. We will 
discuss examples and their topology in Sec. 2. 

The third method consists of taking known solutions and mapping them to new solutions. 
This technique is similar in spirit to conformal mapping. Any analytic function of one 
complex variable, w = f(z) maps z = x + iy to w = u(x,y) + iv(x,y), where, if i/;(x,y) 
is a solution to Laplace's equation, *$f[u(x,y),v(x,y)] is also a solution. This idea exploits 
the invariance of Laplace's equation in two dimensions under the two-dimensional conformal 
group. If one can find groups under which the Einstein equations are invariant (or under 
which some class of Einstein equations for special systems are invariant), one can generate 
new solutions from old. This technique is of wide applicability, and has been exploited 
heavily in general relativity to find new exact solutions. It is by no means specific to 
cosmological models, but it can be applied in cosmological situations, and has, perhaps, 
been underutilized in this field and deserves more attention there. 

In the specific case of stationary axisymmetric spacetimes, early generating methods such 



as the Kerr-Schild ansatz I8l.l9l ll0l .llll]. the complex transformation method 12J, [l3|, [l^ , and 
Hamilton- Jacobi separability [15[ were used to derive new solutions, but all of them are 
contained in the charged Kerr-Taub-NUT class (with cosmological constant). An im- 
portant development for the derivation of exact stationary axisymmetric solutions was made 
by Ernst 0, Q] , who obtained a new representation of the corresponding field equations 
which is independent of the coordinates chosen and, therefore, allows one to investigate the 
symmetries involved and find new solutions. Ernst also proposed two generating techniques 
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19. 20]. Using the compact Ernst formu 



21 



22 



23, 



24 



that were later enlarged and unified by Kinnersley 

lation of the field equations it was also possible to obtain new solutions 
possessing a certain prescribed polynomial form of the Ernst potential. Modern solution 
generating techniques involve Lie groups of transformations or Backlund transformations. 



The first such group was found by Geroch |26l.l27j] and generalized by Cosgrove |28|,|29j]. The 
Geroch group was investigated very intensively and it was found that it possesses subgroups 
that preserve asymptotic flatness 

00,0, HQ- 

Today, there exist three main solu- 
tion generating tecniques: HKX (Hoenselaers, Kinnersley, Xanthopoulos) transformations 
[34I . which are based upon special subgroups of the Geroch group, Backlund transforma- 




tions [36. l37], w hich are applied directly to the Ernst potential, and the inverse scattering 
, 3 ' which is based on a reformulation of the nonlinear field equations as a linear 
eigenvalue problem. Later on it was shown j^oi 41 1 that all these techniques are related to 



one another, and can be used to generate the same type of solutions (for an introductory 
review of the solution generating techniques and the relationships between them see [ij). 

The mapping method is by no means specific to stationary axisymmetric metrics, but 
it can be applied to any spacetime characterized by two or more commuting Killing vector 
fields as in the case of the cosmological models investigated in this work. We will give some 
examples in Sec. 3 and a guide to possible new uses of the technique. 

The paper is organized as follows. Section 2 will give examples of horizon and causal 
structure methods and solutions obtained using them. Section 3 will give one example of 
mapping methods and discuss others. Finally, Section 4 will discuss possible new directions 
in the use of known solutions to find new cosmologies. 



II. HORIZON AND CAUSAL STRUCTURE METHODS 



A. Horizon Methods 



There are a large number of solutions of the Einstein equations that have cosmological 
coordinate patches as well as "black hole" regions that are separated by horizons. The major 
problem with the cosmologies generated in this way is they are incomplete manifolds, and 
what have in the past been interpreted in some of them as singularities are just the horizons 
that separate one part of the complete manifold from the other. A second problem, as we 
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have seen, is that in certain cases the complete manifold is of such high symmetry that is is 
impossible to distinguish between the cosmological region and the "black hole" region. The 
most blatant example of the second problem is just ordinary Minkowski space. 

If we return to Figure ^ we can think of this diagram as a picture of Minkowski space 
with the meeting point of lines aa and bb an arbitrarily chosen origin O. The usual metric 
is 

ds 2 = -dT 2 + dx 2 + dy 2 + dz 2 . (8) 

In region I the hyperbolic lines are Just the usual orbits of a particle with constant accelera- 
tion that form the basis of Rindler j43( space, and if the lines are given by x = \J1r — 1 cosh t 
and T = \J1r — 1 sinht, z = y = const., the metric becomes 

ds 2 = -{2r-l)dt 2 + ^—^ + dy 2 + dz 2 . (9) 

In this coordinate system the metric is very similar in form to the Schwarzschild metic, and 
there is a coordinate singularity at r = 1/2, which is a horizon similar to that of the deSitter 
metric. 

In Figure ^ the line bb is the surface r = 1/2, and for r < 1/2 we can make the t <-> r 
coordinate change and the metric becomes 

ds 2 = -—^-dt 2 + ( X ~ 2t ) dr2 + + d * 2 > ( 10 ) 

a cosmological model. If we change coordinates using t = ~[1 — r 2 {cosh 2 p — 
sinh 2 p sin 2 9 cos 2 </?}], r = tanh -1 [tanh p sin 9 cos ip], y = r sinh p sin 9 sin ip, z = r sinh p cos 9, 
the metric (|1U|) becomes the Milne universe, a flat cosmological model with metric 

ds 2 = -dr 2 + r 2 [dp 2 + sinh 2 p{d9 2 + sin 2 9d<p 2 )}, (11) 

a k = — 1 FRW metric. Of course, the manifold is nothing more than flat space and the 
breaking up of the manifold into regions I-IV is an observer-dependent phenomenon due 
to singling out the origin O as the position of a special observer, while there is actually 
no physical reason that this point is priveleged over any other point of the manifold, and 
there is no real difference between any of the four regions. Basically, we can say that the 
same is true of deSitter space, which is just a constant curvature space everywhere, and the 
"cosmological" region has nothing to distinguish it from any other region of the manifold. 



S 



There are a number of metrics which are not as degenerate as the Minkow- ski and 
deSitter cases whose cosmological regions are physically distinguishable from the "black 
hole" regions, and the resulting cosmological models have long been known and are named. 
Perhaps the simplest of these is the Kantowski-Sachs-Schwarzschild manifold. Here we can 
still use the generic Figure Q to represent the Penrose diagram of this metric, with region I 
representing the Schwarzschild coordinate patch where the hyperbolic lines are r = const, 
lines and the lines bb is the horizon at r = 2m. In region II the hyperbolic lines are still r = 
const, lines, but they no longer represent spacelike surfaces. For r > 2m the metric can be 
written as 



but for r < 2m the Killing vector d/dr is no longer spacelike, and we may make the 
transformation r <-> t mentioned in Sec. 1, and we find 



which is an obvious cosmological model, the vacuum Kantowski-Sachs model. The fact that 
at r = (t = in the new coordinate system) which is represented by the line at +F is a true 
curvature singularity distinguishes region II from region I and makes it impossible to find a 
global coordinate transformation that makes this cosmological model indistinguishable from 
the "black hole" region. Of course, however, this cosmological model, and all cosmological 
models generated by horizon methods, have the problem that what has been regarded as 
singularities in these models may not be a curvature singularity, but only a horizon where 
we pass from one coordinate patch to another. In the case of the Kantowski-Sachs solution 
the t = singularity is one of the usual singularities of the model, and it is a true curvature 
singularity. The other "singularity" at t — 2m is only a null surface. 

While the Kantowski-Sachs solution, even though it was originally discovered by means 
of studies of groups of motion that were not transitive on three surfaces, and constituted 
a generalization of the Bianchi cosmological models, and was immediately recognized as 
the Schwarzschild solution inside the horizon, in the next model we will discuss the two 
regions on either side of the horizon were studied separately. The "black hole" region was 
the NUT space of Newman, Tamburino and Unti £9j, and the cosmological region was the 
Taub cosmology [44J. Misner [4^ showed that these two solutions could be seen to be part 
of a larger manifold (up to topological questions which will be discussed below). In the black 




(12) 




(13) 
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hole region one can write the NUT metric as 

y2 27717" I 2 T 2 -\- I 2 

ds2 = TTT? (^ + 2/cos^) 2 + -dr 2 

r z + l z r z — 2mr — l z 

+ (r 2 + I 2 )(d8 2 + sin 2 9d<p 2 ). (14) 



For the region of this metric where < r < m(l + J 1 + l 2 /m 2 ), the r <-> t transformation 
gives 

j2 l 2 + 2mt-t 2 . J oJ 2 t 2 + l 2 „ 

ds = 5 ( dr + 21 cos Udu) ) —— -dt 

t 2 + l 2 v rj l 2 + 2mt-t 2 

+ (t 2 + l 2 )(d6 2 + sin 2 Odtp 2 ), (15) 

which is an anisotropic cosmological model. The major problem with this model is that 
in the cosmological sector we can rewrite the metric on t = const, surfaces by making the 
coordinate transformation 

V = - tt/2 , (16) 
z = r/2l , (17) 
x = <p, (18) 



which gives the three-dimensional line element as 

+ 2mt 
t 2 + l 2 



Al 2 (l 2 + 2nd - t 2 ) 

da 2 = — - — -(dz - sin ydx) 2 + (t 2 + l 2 )(dy 2 + cos 2 ydx 2 ), (19) 



which can be written as 

do* = ^l 2 +^rnt-t 2 ) {u , ? + (f2 + l2)[{ulf + ^ (2Q) 

where u 3 = dz — sinydx, u 1 = cos y cos zdx — sin zdy, u 2 = cos y sin zdx + cos zdy. These 
differential forms are the invariant one-forms of the Bianchi type IX cosmological models, 
the invariant one-forms on S 3 . This means that the "natural" topology of the three-surfaces 
given by (J2U|) is that of S 3 , with < z < Arc, < x < 2tt, < y < n. If one maintains 
this topology on passing through the horizon, the vector d/dz becomes timelike, and the S 3 
topology implies the possibility of closed timelike lines. This topological difficulty is present 
in many cosmological models generated by horizon methods, as we will see below. 

The metrics we have discussed seem to be the only ones where the cosmological models 
have been noticed explicitly to be a part of a larger manifold that has non-cosmological 

10 



sectors. For example, it is only recently that the Kerr metric inside its horizon has been 
considered as a cosmological model. Of course, the metric in this part of the Kerr manifold 



has been studied 



471 ] , but there seems to have been no attempt to identify it as a cosmological 



model. In Boyer-Lindquist coordinates the Kerr metric has the form 

,9 r 2 — 2Mr + a 2 . , . 9 „ , ,, 9 sin 2 9 r , 9 9 . , , , l9 

ds = 9 9 T7T i dt -asm 2 Sty] 2 + — - — r 2 + a 2 )d<j) - adt} 2 

r z + a z cos z 9 r z + a z cos z 9 

+ r 2 + fj°^ 9 2 dr2 + ( r2 + a " cos2 6 ) d ° 2 ( 21 ) 
r z — 2Mr + a z 



This metric has two horizons, r± = M±\ / M 2 — a 2 , and beyond the outer horizon, r + , d/dr 
is a spacelike vector, and the metric represents a spinning black hole. At r + the light cones 
tip over to the point where d/dr becomes timelike and we can make the transformation 
t <-> r. Unfortunately, at the inner horizon r_ the light cones tip back to the point where 
d/dr becomes spacelike again. In the region r_ < r < r + the Kerr metric becomes a 
cosmological model, 

, 9 2Mt — t 2 — a 2 t1 . n „ , ,,n sin 2 6 r . 9 9 . , , , l9 

ds = ~T9 9 T7T i dr - asm 2 0# 2 + — t 2 + a 2 )d(j) - adr} 2 

V + a z cos z 6 t z + a z cos z V 

f2 + " 2co * 29 dt 2 + (t 2 + a 2 cos 2 9)d6 2 . (22) 



2Mt -t 2 -a 2 
This metric, with the simple transformation [3], 



t = a yi- /3 2 cos{e- T ) + l], (23) 
a = M, (3 = a/M (0 < f3 <l) transforms into 



ds 2 = e- x/2 e T/2 (-e- 2T dr 2 + d9 2 ) + ayjl - (3 2 sin( e - T ) [e p d6 2 

+2e p Qd5d<j) + {e p Q 2 + e~ p sin 2 6)d<j) 2 ], (24) 



where 



A = t - 2 ln(a 2 { [J I - f3 2 cos{e~ T ) + l] 2 + (3 2 cos 2 6}) , (25) 



P = ln[(l - (3 2 ) sin 2 (e- T ) + f3 2 sin 2 9] - \n[a^l - (3 2 sm( e - T )} 

- ln([ v /l-/3 2 cos(e- T ) + l] 2 + (3 2 cos 2 9), (26) 

2a/?sin 2 g[v^^cos(e-) + l] 
^ (l-/5 2 )sin 2 (e-)+/3 2 sin 2 ' J 

The metric (J2U) is a Gowdy cosmological model with S* 1 x S 2 topology 
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FIG. 2: The evolution in r of P — h as a function of 9 from Eqs. (|25II27[) for a = 1, /? = 1/2. 
Fig. (a) corresponds to r = — ln(7r), (b) to r = —1.144, (c) to r = —1.1, (d) to r = —0.75, (e) to 
r = +5, (f) to r = +10 

The solution given in Eqs. (|231 -|2Zj) is a function of 9 for < 9 < n and is valid for 
values of r that correspond to the part of the manifold between r_ and r + , that is, for 



— ln(7r) < t < +oo, (note that t takes the values M ± a/ M 2 — a 2 at the two limiting values 
of r). 

As an example, Figs. El and H3 give P — h (h = — lnfa^l — sin(e _T )]) and Q as functions 
of ^ for various values of r for a — 1, /3 = 1/2. The evolution of these two functions is that 
of a "spike" that at the limiting values of r is practically flat (except at 9 = 0, tt, where it 
drops off drastically), and which becomes sharper for intermediate values of r. It should 
be mentioned that the Kerr metric is a special case of what is usually called the Kerr-Taub- 
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FIG. 3: The evolution of Q as a function of 9 in r. Fig. (a) corresponds to r = —1.1439, (b) to 
r = -1.1, (c) to r = - ln(vr/2), (d) to r = +5 

NUT metric, which is geven in Ref. H- 

This metric has three parameters, a, m, I, the 
Kerr parameter, the mass, and the NUT parameter respectively. For a = we have the 
metric 1)140. and for I = we have the Kerr metric. This means that the Gowdy model 
we have given is also a special case of the "Taub" region of Kerr-Taub-NUT. However, the 
Taub region of this manifold, in spite of its name, never seems to have been considered as a 
cosmological model. 

As we have stated several times, the models studied so far are vacuum models. While 
matter-filled models may, in principle, be generated by horizon methods, in at least some 
cases the models with matter have very different behavior from the vacuum case, and may 
preclude naive generalizations of the method to the non-vacuum case. For instance, it is 
known that the Taub model with fluid matter has true curvature singularities where the 
vacuum case has null surfaces |49l. The same is true for the Kantowski-Sachs model 1501. 
However, the Brill model [5J| is a Taub-NUT electrovac solution which is non-singular in 
the same way as the Taub-NUT vacuum case. In fact, the Reissner-Nordstr0m solution with 
electric charge large enough has two horizons, and between the outer and inner horizons 
there is a region that can be interpreted as an electrovac cosmology that has no curvature 
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or electromagnetic field singularity at the horizons. 

Besides the solutions we have mentioned, there are numerous metrics that represent 
gravitational fields outside some material body which in vacuum have "black hole" sectors 
outside of some horizon. Of course, no hair theorems [52] say that there are no true vacuum 
black holes that are not Kerr or Schwarzschild. Solutions of the type mentioned above seem 
always to have curvature singularities on or outside their horizons where they can be seen by 
observers at infinity. While this fact makes them poor candidates for black holes, it makes 
them interesting as cosmological models inside their horizons. The "singularities" of the 
models with curvature singularities on their horizons will be inhomogeneous, with part of 
the singularities being simple horizons, but other regions will be true curvature singularities. 
These structures might be of great interest in the study of cosmological sing ularities. Perhaps 
the simplest of the metrics of this type might be the Tomimatsu-Sato 0^ 13 models which 
do have inhomogenous horizons, but there are other known solutions, for example, one with 
many multipoles and large curvature singularity regions on the horizon js^ . 



B. Causal Structure Methods 



This method seems to have been little used, and perhaps the most telling reason for this 
is that there are very few examples where it works. The only example that has wide currency 
is that of certain Gowdy models. If one begins with the Einstein- Rosen waves j^, which 
have the metric 

ds 2 = e 2 ^(dr 2 - dt 2 ) + rV 2 W + e^dz 2 , (28) 
with 7 and ip functions of r and t, the Einstein equations for this metric are 

+ -VV - Att = 0, (29) 
r 

%r = r(ij>*+ij>% 7, t = 2rWV> t , (30) 

where it is well known that the equations for 7 can be integrated directly once ip is known, 
since the integrability condition for these two first-order partial differential equations is just 

If we now simply make the coordinate change t <-> r, we have 

ds 2 = e 2 ^(dt 2 - dr 2 ) + t 2 e- 2 ^d<p 2 + e 2 ^dz 2 . (31) 
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Here the only obstacle to interpreting this metric as a cosmological model is the wrong signs 
in the dr 2 and dt 2 terms if we want to interpret d/dr and d/dt as spacelike and timelike 
vectors respectively. If we add a complex constant to 7, 7 = 7 + iir, we change the signs 
of g tt and g rr , and since we have done nothing more than change the names of r and t, the 
Einstein equations are unchanged (since they depend only on derivatives of 7), and we have 



These are the original equations for the unpolarized Gowdy T 3 cosmological model p . 

Notice that in order to change the sign of the dt 2 — dr 2 term, it was necessary to add 
— 171 to 7, and if the Einstein equations were to depend on 7 in any other way than through 
derivatives or e 7 , they would become complex and then would not necessarily have real 
solutions for 7 and ip. This difficulty is a paradigm for the problems one would encounter 
in trying to use at *-> r type of coordinate change. One would not expect this technique 
to be successful in the great majority of cases, especially if there were matter present. This 
seems to be reflected in the fact that there exist few examples of this method. Note that the 
examples we have given, the plane wave of the Introduction and the T 3 Gowdy model, are 
vacuum spacetimes, where the temporal variable and one of the space variables appear in 
the form of —dt 2 + dz 2 (in the Gowdy model, multiplied by a conformal factor). While we 
will not try to explore this kind of idea further here, perhaps theorems about the existence 
of solutions found by means of the method could be based on this fact. 

C. Topological Questions 

An interesting feature of cosmological models generated by means of horizon and causal 
structure methods is that for philosophical reasons one often wants these models to have 
compact t = const, sections, that is, they should be closed cosmological models. This can 
be achieved for many of the solutions given above by simply specifying the global topology 
without changing the local geometry. This has caused some problems with these models, 
especially the vacuum models generated by horizon methods, since they represent incomplete 
manifolds, and when one passes through the horizon, the spacelike direction in which the 



V 

7, i = t(^ + ^ 2 r ), nr, r = 2tlp,t1p,r 



(32) 



(33) 
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manifold would have to be closed becomes timelike, leading to the possibility of closed 
timelike lines in the "black hole" sector. 

Perhaps the simplest model where this occurs is the Kantowski-Sachs- Schwarzschild 
manifold. If we use (II 3|) for the Kantowski-Sachs model, notice that the 6<p sector has a 
natural two-sphere topology, and there is no obstruction to compactifying in the r direction 
by simply assuming that < r < 2tc, with r = and r = 2ir identified, giving the 
manifold an S 1 x S 2 topology. Of course, if we consider the other side of the horizon, 
the t <-> r reparametrization means that the r-direction becomes the t-direction, and we 
have < t < 2tt and the possibility of closed timelike lines. In principle one could have 
timelike lines which pass from the cosmological region to the black hole region and remain 
trapped in that region in an eternal closed timelike curve. In the case of the Kantowski- 
Sachs- Schwarzschild manifold it seems to be impossible for this to happen. From Fig. 0] 
one can see that timelike geodesies which begin inside the horizon will always stay within 
the horizon, simply passing throu gh t he crossing point of the past and future horizons (a 
focussing point for all geodesies) [55]. The choice of topology of this metric has a long 
history. The original studies of this metric assumed that r ran from — oo to +oo and that 
the topology was R 1 x S 2 . This topology was usually used until Laflamme and Shellard 



introduced the S 1 x S 2 topology 5t\. Since that time most authors have used it, but some 
still prefer the R 1 x S 2 topology 57 1. 

Probably the first solution where topological problems were noticed was the Taub-NUT 
model, where the natural topology of t = const, slices in the Taub model is S 3 , which means 
that the same topology must apply to the tOtp sector of NUT space. Here it is possible for 
timelike geodesies to leave the Taub region and pass into the NUT reg i on, where they may 
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be trapped in trajectories that return to the same spacetime point |45l . 

For the other solutions generated by horizon methods mentioned so far, it seems that we 
will have similar problems, but there has been no study of the effect of topology changes 
on the causal properties of the solutions. Even degenerate solutions such as the deSitter 
model and Minkowski space can be given different topologies by means of this type of 
identifications. Thus, for example, the Milne model, while it is locally flat, may be made 
into a more acceptable cosmology by means of such a global topology change, but it would 
suffer from the same problem as the other metrics of this type in that the Rindler space 
sector would then have an unacceptable topology. The Kerr-Gowdy solution has this same 
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FIG. 4: A Kruskal diagram for geodesies in Kantowski-Sachs-Schwarzschild. For an open topology 
the timelike geodesic F can leave the cosmological region, enter the Schwarzschild region and return 
to the cosmological region. This behavior is similar to that of some geodesies in Taub-NUT. For 
the closed topology, however, geodesies wrap around between r = and r = 2ir (r being the 
new "radial" direction inside the horizon), following the solid-dashed sawtooth path (here, as an 
example, for a lightlike geodesic), which always stays inside the horizon except at the crossing 
point of the past and future horizons 

topological difficulty, since the interior is supposed to have an S 1 x S 2 topology, and outside 
the horizon this implies an S 1 topology in the timelike direction. The causal properties of 
this solution with this topology seem never to have been studied. 

For metrics generated by causal structure methods, there is, in principle no obtruction 
to topology changes, since there is no "black hole" sector where new topologies can lead 
to causality problem. The Gowdy T 3 model is a good example of this. While we have not 
specified the ranges of r, and z in (|3ip. Gowdy, wishing to have a closed topology for the 
t = const, surfaces in his model, assumed that < r, <f), z < 2n with points at and 27T 
identified, which gives the model a three-torus topology. Since there is no "black hole" sector 
for the new manifold, there is no reason not to choose this topology. The only difference 
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that this change introduces is in the boundary conditions that ip and 7 must satisfy. If r, 
<f), and z run from —00 to +00 the solutions of Eq. are built up of eigenfunctions of 

continuous eigenvalues, while in the Gowdy topology the eigenvalues are discrete. 

Note that if we were to try to make the t *->■ r change for this topology, the Einstein- 
Rosen wave thus generated would be closed in the time direction, but the Gowdy manifold 
is complete, and there is no reason to concern ourselves with this possibility. 

As we mentioned in the Introduction, the somewhat artificial model we gave there may 
be closed in the space directions without changing the local metric. 

III. MAPPING METHODS 

As we mentioned in the Introduction, the solution generating tecniques developed in the 
last few decades have not often been used in the context of cosmological models. However, 
it is known that these techniques can be applied when the spacetime possesses at least two 
commuting Killing vector fields. This is exactly the case of the cosmological models under 
investigation in this paper. 

In the original works on solution generating techniques the central idea was to project 
the spacetime on the two-dimensional hypersurface defined by the Killing vector fields. 
This projection allows one to reformulate the spacetime metric and the corresponding field 
equations in such a way that the symmetries of the differential equations can be investigated 
in a straightforward manner. An alternative but related procedure is to derive the Ernst 
representation of the field equations where it is possible to apply (with some modifications) 
known techniques. In the present work we will use this second alternative. 

A. Ernst Representation of the T 3 Gowdy Models 

In this section we will begin by considering the unpolarized T 3 Gowdy model. We will 
use a slightly different parametrization of these models from that used in (|3T|) in order to 
make comparisons with Refs. j^. The unpolarized model has the metric 

ds 2 = e~ x/2 e T/2 (-e- 2r dT 2 + d6 2 ) + e - T [e p (da + Qd5) 2 + e - p d5 2 } , (34) 

where the functions A, P and Q depend on the coordinates r and 9 only, and < a, 5, 9 < 2n. 
These spacetimes are characterized by the existence of two commuting Killing vector fields 
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r/i = 8/ da and rj 2 = 8/ 85. In the special case Q = 0, the fields rji and rj 2 become hypersurface 
orthogonal to each other and the metric (}3"4*|) describes the polarized T 3 Gowdy models. 

Einstein's vacuum field equations for the T 3 models consist of a set of two second order 
differential equations for P and Q 

P TT - e- 2 ^P fi e ~ e 2P (Q 2 T - e^Q 2 ,) = , (35) 

q tt _ e - 2T Q, ee + 2(P T Q, T - e' 2r P fi Q,e) = , (36) 
and two first order differential equations for A, 

A)T = P 2 + e~ 2T + e 2P (Q 2 T + e^Q 2 ,) , (37) 



A, fl = 2(P fl P T + e 2T Q,,Q iT ) . (38) 

The set of equations for A is the equivalent of Eqs. (jHB^) and can be solved by quadratures 
once P and Q are known. 

As we mentioned, the special case of the polarized T 3 model is obtained from the metric 
(JHU) just by taking Q = 0. The resulting field equations are easier to handle, and a general 
solution for the main function P can be found by separation of variables. In fact, let us 
consider P(t,9) = T(r)O(0); then, Eq. (|35j) separates into 

ld 2 Q „ , d 2 T 



—n 



and — + n 2 e~ 2T T = , (39) 



6 d6 2 ' ~ dr 2 

where n is the separation constant which in this case has to be an integer in order for the 
condition 0(0 + 2ir) = 0(0) to be satisfied. With this assumption, the general solution for 
P can be written as an infinite series of the form 
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P = ]T [A n cos(n0) + B n sin(nfl)] [C n J (n e - T ) + D n N (ne~ T )] , (40) 

n=0 

where A n , B n , C n and D n are arbitrary constants. If we want to avoid singularities at 
r = +oo, the constant D n has to vanish. 

The models described by the general (Q ^ 0) metric (|34j) have been used extensively for 
numerical investigations in classical as well as in minisuperspace quantum gravity |k| . One 
of the reasons why these investigations have used numerical methods is because it is usually 
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believed that the set of main field equations (|3*5Jl and (jHUj) is such a complicated system that 
analytic solutions would be difficult to find. We will show here that it is possible to generate 
unpolarized solutions (Q ^ 0) from a given polarized solution {Q = 0) by using modern 
solution generating techniques. Of course, the new solutions will be particular solutions, 



and the numerical investigations in Refs. 



59| were carried out with the aim of studying the 



general behavior of the models near a singularity, information that no particular solution 
can give. However, families of particular solutions can give us clues about how to set up 
numerical solutions. 

We can apply the solution generating techniques by writing the field equations in such 
a way that the symmetries involved can be derived and understood easily. To this end, we 
first introduce a new "time" coordinate t = e~ T and a new function R = R(t, 6) by means 
of the equations 

R,t = te 2P Q fi , Re = te 2P Q, t . (41) 
Then, the field equations (}3T)j) and (J33j) can be expressed as 

t 2 (P,u + \p,t ~ P fi e) + e~ 2P (Rl - R 2 e ) = , (42) 

te p (R, tt + jR, t - R fi e) ~ 2[(te p ), t i? it - (te p ) fi R,e] = . (43) 

Furthermore, this last equation for R turns out to be identically satisfied if the integrability 
condition R te = R et is fulfilled. 

We can now introduce the complex Ernst potential e and the complex gradient operator 
D as 

e = te p + tR J and D = ^ , z^J , (44) 
which allow us to write the main field equations in the Ernst-like representation 

Re{e) (D 2 e + l -Dt De^j - {Def = . (45) 

It is easy to see that the field equations (j42[ I4*3|) can be obtained as the real and imaginary 
part of the Ernst equation (|43|) . respectively. 

The importance of this representation is that it also can be derived from a Lagrangian 
by means of a variational principle. In turn, that Lagrangian may be interpreted as a 
metric Lagrangian defined in a two dimensional Riemannian space (the potential space), 
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the coordinates of which are the real and imaginary parts of the Ernst potential. Applying 
the variational principle in the potential space on the metric Lagrangian, one obtains the 
corresponding geodesic equations, which turn out to be equivalent to the Ernst equation 
(J45|) . Hence, to investigate the symmetries of the Ernst equation one can study infinitesimal 
transformations which leave invariant the geodesic equations in the the potential space. In 
particular, the transformations associated with the Killing vector fields of the metric in the 
potential space leave the corresponding geodesic equations invariant. One could think of a 
solution of the Ernst equation as a geodesic in the potential space, and the Killing vectors of 
the metric as transformations that starting from a given geodesic lead to a different geodesic, 
i.e., to a different solution of the Ernst equation. This is the basic idea behind some of the 
known solution generating techniques. 

We will now derive a simple but illustrative symmetry of the Ernst equations which allows 
to generate new solutions. To this end, we introduce a new complex potential £ = 6) by 
means of the relationship 

«-£f (46) 
Then, the Ernst equation ()45|) transforms into 

(1 - iC) (p*£ + \DtDdj + 2£*(£>£) 2 = , (47) 

where represents the complex conjugate potential. Furthermore, we introduce new coor- 
dinates x and y by 

t 2 = c 2 (l - x 2 )(l - y 2 ) , 6 = cxy , (48) 

where c is a real constant. In these coordinates, the main field equation (}4Tj) can be written 
in the following form 

(1 " ££*){[(! - ~ K 1 - y 2 )tvlv} + 2 C[(1 - x 2 )i 2 x - (1 - y 2 )i 2 y ] = , (49) 

which explicity shows the invariance with respect to the change of coordinates x <-> y, i.e., 
if £,(x,y) is a solution of (|4*D^) then £(y,x) is also a solution. The simplest solution of Eq. 
(j4Tlj) is = x, so £ _1 = y is also a solution. A linear combination of these two solutions 

= ax + iby (50) 

turns out to also be a solution if the condition a 2 + b 2 = 1 is satisfied. In this way it is 
relatively easy to generate new solutions. Once the potential £ is known one can calculate 

21 



the functions P and R algebraically by means of Equations (|46J) and (|44jl . Then, Eq. ()41|) 
may be integrated to yield Q. Finally, the solution for the function A may be calculated by 
quadratures. 

It is interesting to note that one can derive some of the main properties of the Gowdy 
model just by inspecting the explicit form of the Ernst potential. If the Ernst potential 
e is real, then the function R vanishes and, therefore, the metric function Q is a constant 
which can be absorbed by means of a suitable coordinate transformation. Thus, a real Ernst 
potential corresponds to a polarized Gowdy model. In the case of complex Ernst potential 
with nontrivial imaginary part (not proportional to its real part), it is guaranteed that the 
resulting metric corresponds to an unpolarized Gowdy model (Q ^ 0). For instance, even 
the simple solution (|5U|) will lead to a nontrivial unpolarized solution with Q ^ 0. 

B. Generation of New Solutions 

In this section we will briefly describe the way new solutions can be generated by using 
HKX transformations. Let us assume that a solution for the polarized Gowdy model is given 
by P = P (t,9) and A = Xo(t, 9). The corresponding Ernst potential is real, and we denote 
it by e = te p °. 

A HKX transformation acting on eo generates a new complex Ernst potential e which will 
correspond to an unpolarized Gowdy model only if its imaginary part is not proportional 
to its real part. This will depend on the explicit form of the initial Ernst potential eo- 
However, if we apply two different HKX transformations it can be shown that the resulting 
Ernst potential is nontrivial. For this reason, we present here the result of the action of two 
HKX transformations on the real Ernst potential eo- For the sake of simplicity, we use the 
coordinates x and y defined in Eq. (J48j) in which the new Ernst potential can be written as 

Q 

_ x(l - hi(jl 2 ) + vyjpa + Mg) + (1 + /fi/fg) + - Mg) /r-iN 
°x(l - fafjv) + iy(n\ + yu 2 ) + (1 + HifMi) - - fa) ' 
where we have introduced new functions \i\ and \i 2 defined by 

/it = aie 2/3 ~ and /i 2 = a 2 e 2l3+ . (52) 

Here oti and a 2 are real constants and j3± are functions which satisfy a set of two first order 
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differential equations 



(x T y){P±)* = (1 T ^)— T (1 - y 2 )^ , (53) 

eo e 



(x T yW±), y = (1T a*/)- 2 * ± (1 - x 2 )^ . (54) 

eo e 

Thus, the generation of new solutions reduces to the integration of the differential equa- 
tions for f3±. The new constants a.\ and a.2 have been introduced by the two HKX transfor- 
mations. For vanishing ol\ and 02, the potential (JoTj) reduces to the original Ernst potential 
of the polarized Gowdy model. 

Equations (|53~j) and (p>fj) allow us to generate new unpolarized solutions starting from a 
given polarized Gowdy model. The derivation of explicit solutions implies the integration 
of a set of two first order differential equations for a function which determines the Ernst 
potential of the new solution. In turn, from the Ernst potential one can obtain the new 
metric functions which completely determine the spacetime of the new cosmological model. 
In order to complete this procedure it is necessary to carry out lenghtly but straightforward 
calculations. We are attempting to calculate the unpolarized solution corresponding to the 
general solution ()40j) of the polarized T 3 Gowdy model. 



IV. CONCLUSIONS AND SUGGESTIONS FOR FURTHER RESEARCH 

We have given several examples of the three methods listed in the Introduction. It is 
obvious that these three approaches are far from exhausted, and may give a number of 
interesting exact solutions to known types of cosmological models and suggest new types 
of models of interest. It is probable that the horizon and mapping methods may be more 
productive than causal structure methods. 

It seems that horizon methods have barely touched the surface of possible solutions. 
There are a large number of exact solutions of "black hole" type given in Ref. 4|| that 
are candidates for generating cosmological models inside their horizons. We have looked 
at the Kantowski-Sachs-Schwarzschild, Taub-NUT and Kerr manifolds in some detail. We 
have also mentioned as possibilities the Tomimatsu-Sato and Quevedo class of solutions 
which should lead to models with mixed singularities, part null surface and part curvature 
singularity. There are large classes of vacuum type D solutions that might give interesting 
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cosmologies [46j. Even such well known solutions as Kerr-Taub-NUT do not seem to have 
been investigated thoroughly as cosmological models. 

Models with matter seem hardly to have been touched. While fluid models may be 
difficult to generate, the example of electrovac universes shows us that there are many 
possible cosmologies of this type. Even the interior of the Reissner-Nordstr0m solution 
does not seem to have been investigated in this context. It would also be interesting to 
investigate the possibility of obtaining useful cosmological models from the most general 
type D electrovac solutions with cosmological constant, found by Debever et al. [6l| . which 
contains 13 different parameters and includes the Kerr- Newman solution as a special case. 
Even the Kerr-Newman solution would generate an R 1 x S 2 (and even, perhaps, an 
S 1 x S 2 ) Gowdy solution, equivalent to, or a generalization of, the solutions of Carmeli, 
Charach and Malin [60] for T 3 models. Other types of solutions, such as scalar field models 
do not seem to have been investigated at all in this context. 

If anything, mapping methods seem even more underutilized in cosmology. The example 
we have seen is one of many that come to mind. For the Gowdy models one could use 
these methods to investigate the S 1 x S 2 models and perhaps generate unkown solutions. 
The corresponding Ernst potential of the field equations present some technical difficulties 
related to the specific topology of the model. This problem is currently under investigation. 
For this case, it would also be interesting to generate new unpolarized solutions starting 
from the general polarized solution, which in the case of S 1 x S 2 models, just as in the T 3 
models, can also be represented as an infinite series of eigenf unctions since the corresponding 
field equations reduce to separable linear second order differential equations. 

Solutions with with true curvature singularities on possible inner and outer horizons 
might be especially interesting. An example of this idea was given by Moncrief (3], who 
used the Geroch group to generate a solution that had a curvature singularity in place of a 
horizon in the Kerr-Taub-NUT case. 

There will certainly be many opportunities to apply all three of the techniques discussed 
in the article to the generation of new cosmologies. 



24 



Acknowledgments 

We wish to thank T. Jacobson, V. Moncrief and 0. Obregon for useful discussions. This 
work was supported in part by CONACyT grant 3567E, and DGAPA-UNAM grants 121298 
and IN106097. 



[1] Krasinski A. (1997) Inhomogeneous Cosmological Models. Cambridge U. P., Cambridge. 

[2] The literature on the deSitter metric is vast. For the description given here see: Hawking 

S., Ellis G. (1973) The Large Scale Structure of Space-Time. Cambridge U. P., Cambridge; 

M0ller C. (1952) The General Theory of Relativity. Oxford U. P., Oxford. 
[3] deSitter W. (1917) Mon. Not. Roy. Astron. Soc. 78, 3. 
[4] Lemaitre G. (1925) J. Math, and Phys. (M.I.T.) 4, 188. 
[5] Robertson H. (1928) Phil. Mag. 5, 835. 

[6] Gowdy R. (1971) Phys. Rev. Lett. 27, 826; (1974) Ann. Phys. (N.Y.) 83, 203. 

[7] Bondi H., Pirani F., Robinson I. (1959) Proc. Roy. Soc. London A251, 519; Ehlers J., Kundt 

W., (1962) in: Witten L. (Ed.) Gravitation: An Introduction to Current Research. Wiley, 

New York. 

[8] Kerr R. P. (1963) Phys. Rev. Lett. 11, 237. 

[9] Newman E. T., Tamburino L. and Unti T. (1963) J. Math. Phys. 4, 915. 
[10] Kinnersley W. (1969) J. Math. Phys. 10, 1195. 
[11] Plebahski J. F. and Demiahski M. (1976) Ann. Phys.(N. Y.) 98, 98. 
[12] Newman E. T. and Janis A. I. (1965) J. Math. Phys. 6, 915. 
[13] Newman E. T. et al. (1965) J. Math. Phys. 6, 918. 

[14] Demiahski M. and Newman E. T. (1966) Bull. Acad. Pol. Sci., Ser. Math. Astron. Phys. 14, 
653. 

[15] Carter B. (1968) Comm. Math. Phys. 10, 280. 

[16] Demiahski M. (1973) Acta Astronom. Pol. 23, 197. 

[17] Ernst F. J. (1968) Phys. Rev. 167, 1175. 

[18] Ernst F. J. (1968) Phys. Rev. 168, 1415. 

[19] Kinnersley W. (1973) J. Math.. Phys. 14, 651. 



25 



[20] Kinnersley W. (1977) J. Math. Phys. 18, 1529. 

[21] Tomimatsu A. and Sato H. (1972) Phys. Rev. Lett. 29, 1344. 

[22] Tomimatsu A. and Sato H. (1973) Prog. Theor. Phys. 50, 95. 

[23] Yamazaki M. and Hori S. (1977) Prog. Theor. Phys. (Lett.) 57, L696. 

[24] Yamazaki M. (1977) Prog. Theor. Phys. 57, 1951. 

[25] Yamazaki M. (1977) J. Math. Phys. 18, 2502. 

[26] Geroch R. (1971) J. Math. Phys. 12, 918. 

[27] Geroch R. (1972) J. Math. Phys. 13, 394. 

[28] Cosgrove C. M. (1977) J. Phys. A: Math. Gen. 10, 1481. 

[29] Cosgrove C. M. (1977) J. Phys. A: Math. Gen. 10, 2093. 

[30] Kinnersley W. and Chitre D. M. (1977) J. Math. Phys. 18, 1538. 

[31] Kinnersley W. and Chitre D. M. (1978) J. Math. Phys. 19, 1926. 

[32] Kinnersley W. and Chitre D. M. (1978) J. Math. Phys. 19, 2037. 

[33] Kinnersley W. and Chitre D. M. (1977) Phys. Rev. Lett. 40, 1608. 

[34] Hoenselaers C, Kinnersley W., and Xanthopoulos B. C. (1979) Phys. Rev. Lett. 42, 481. 

[35] Hoenselaers C, Kinnersley W., and Xanthopoulos B. C. (1979) J. Math. Phys. 20, 2530. 

[36] Harrison B. K. (1978) Phys. Rev. Lett. 41, 1197. 

[37] Neugebauer G. (1980) J. Phys. A: Math. Gen. 13, 1737, L19. 

[38] Belinsky V. A. and Zakharov V. E. (1978) Soviet Phys.-JETP 48, 985. 

[39] Belinsky V. A. and Zakharov V. E. (1978) Soviet Phys.-JETP 50, 1. 

[40] Cosgrove C. M. (1980) J. Math. Phys. 21, 2417. 

[41] Cosgrove C. M. (1982) in: Proceedings of the Second Marcel Grossmann Meeting on General 

Relativity, Ruffini R. (Ed.), North-Holland, Amsterdam. 

[42] Quevedo H. (1990) Fortschr. Phys. 38, 10. 

[43] See, for example, Rindler W. (1969) Essential Relativity: Special, General, and Cosmological 

(Second Ed.). Springer- Verlag, Heidelberg. 

[44] Taub A. (1951) Ann. Math. 53, 472. 

[45] Misner C. (1963) in: Relativity Theory and Astrophysics 1. Ehlers J. (Ed.), American Math- 
ematical Society, Providence R. I.; Misner C, Taub A. (1968) Zh. Eksp. Teor. Fiz. 55, 233 
(Soviet Physics - JETP 28, 122 [1969]). 

[46] Kramer D., Stephani D., Herlt E., MacCallum M., Schmutzer E. (1980) Exact Solutions of 



26 



Einstein's Field Equations. Cambridge U. P., Cambridge. 
[47] See Burko L., Ori A. (Eds.) (1997) Proceedings of the Workshop on the Internal Structure of 

Black Holes and Spacetime Singularities. IOP, Bristol. 
[48] Obregon O., Ryan M. (1998) IgTqc/ 981 0068 1 

[49] Ryan M., Shepley L. (1975) Homogeneous Relativistic Cosmologies. Princeton U. P., Prince- 
ton. 

[50] Conradi H-D. (1995) Class. Q. Grav. 12, 2423. 
[51] Brill D. (1964) Phys. Rev. 133, B845. 

[52] Israel W. (1967) Phys. Rev. 164, 1776; Carter B. (1970) Phys. Rev. Lett. 26, 331; Hawking 

S. (1972) Commun. Math. Phys. 25, 152. 
[53] Quevedo H. (1986) Phys. Rev. D 33, 324. 
[54] Einstein A., Rosen N. (1937) J. Franklin Inst. 223, 43. 
[55] Jacobson T. (1998) Private communication. 
[56] Laflamme R., Shellard E. P. S. (1987) Phys. Rev. D 35, 2315. 

[57] Xanthopolous B., Zannias T.(1992) J. Math. Phys. 33, 1415; (1991) J. Math. Phys. 30, 1121. 
[58] Hawking S., Ellis G. in Ref. |2|. 

[59] See Berger B., Moncrief V. (1993) Phys. Rev. D 48, 4676; Berger B., Garfinkle D. (1998) 

Phys. Rev. D 57, 4767, and references therein. 
[60] Charach Ch. (1979) Phys. Rev. D 19, 1058; Phys. Rev. D 19, 3516; Carmeli M., Charach Ch., 

Malin S. (1981) Phys. Reports 76, 1. 
[61] Debever R., Kamran N., McLenaghan R. G. (1984) J. Math. Phys. 25, 1955. 
[62] Moncrief V. (1987) Class. Q. Grav. 6, 1555. 



27 



